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1 Introduction 

In this contribution we review a method for deriving collective modes of pair-correlated 
neutron matter as found in the inner crust of a neutron star. We discuss two classes of 
modes lying either above or below the pair-breaking continuum limit 2A and having 
energy spectra of the form e = ujq + aq^, where q is the wave vector of the mode 
and cuo = 2A is a threshold frequency determined by the pairing gap A. One class of 
modes, corresponding to density oscillations, arises in the spectrum for small q when 
a < 0. These modes can be associated with undamped particle-hole bound pairs - ex- 
citons - existing in superfluid neutron matter. We also discuss the modes with a > 0. 
These are spin diffusive modes that arise in neutron-star crusts due to spin perturba¬ 
tions and are located above the pair-breaking threshold. In contrast to excitons, these 
modes are damped. As an application we compute the thermal conductivity due to 
the spin diffusive modes existing in the crust of a neutron star. The general formalism 
that we introduce below is relevant for a number of fermionic ensembles experiencing 
attractive interactions, which become superfluid or superconducting below a critical 
temperature Tc. Examples are ultra-cold atomic vapors, nuclear and neutron matter, 
and deconhned quark matter (for review articles see Ref. [1]). 

Microscopic description of the low-energy dynamics of superfluid systems can be 
formulated within the propagator formalism, where the central quantity is the re¬ 
sponse function. Several types of response functions may be considered, depending 
on the problem of interest. For neutral fermionic systems, the central quantities are 
the density and spin responses. For charged superconductors the charge current re¬ 
sponse has been especially important historically, since it describes the experimentally 
observed Meissner effect. In the present work we focus instead on the pair-breaking 
response of a neutral, non-relativistic superfluid exhibiting A-wave pairing, a response 
which is dominant at temperatures slightly below the critical temperature of the su¬ 
perfluid phase transition. 
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Collective modes of neutron-star inner crusts, which are composed of an ionic 
lattice, a charge-neutralizing electron background, and a superfluid neutron fluid, have 
been studied from a number of perspectives. The electron-ion system oscillates at the 
plasma frequency in the absence of the neutron fluid [2]. However, coupling of the 
lattice phonons to the neutron fluid renormalizes the spectrum of lattice oscillations 
[3] and may lead to instabilities |1]. The low-lying oscillatory modes with spectra of 
acoustic type, notably the Anderson-Bogolyubov mode and the lattice phonon mode, 
tend to mix since they cover the same energy range. Their properties and mixing have 
been investigated in Refs. [HI El [TJ [HI E] . In a first approximation, the pair-breaking 
modes can be treated separately, as they arise near the pair-breaking threshold 2A. 
Before discussing concrete applications to neutron-star crusts, we first outline the 
general formalism for deriving these modes, following closely Ref. [8]. 


2 Formalism 

The techniques we employ were first developed by Abrikosov and Gor’kov (AG) to 
study the Meissner effect in the electrodynamics of dirty superconductors ra- In this 
theory the electromagnetic response of the system is expressed in terms of Matsub- 
ara imaginary-time propagators, assuming a contact pairing interaction. From the 
formal standpoint, the AG theory can be viewed as a reformulation of the theories 
of collective excitations in superconductors developed by Bogolyubov, Anderson, and 
subsequently others using the second-quantized operator formalism. The ideas of 
Landau Fermi-liquid theory (which is applicable to fermionic systems with arbitrarily 
strong interactions) were extended to superfluid systems by Larkin and Migdal and 
by Leggett mm- They reproduced the basic results known at the time, going on 
to derive the so-called Fermi-liquid corrections that arise in strong-coupling theories. 
In contrast to the AG treatment, the latter approaches distinguish particle-hole (ph) 
and particle-particle (pp) interactions having different strength/or sign, in the spirit 
of the Fermi-liquid theory. There is renewed interest in studying the response of su¬ 
perfluid/superconducting nuclear matter, in part with the objective of understanding 
the collective modes of such systems (see Ref. [8] and works cited therein). 

In contributing to this trend, we begin by considering R-wave paired neutron 
matter and write the interactions in the particle-particle and particle-hole channel in 
the forms m 

fap-yS - /p^(*Cr2)a/3(fa2)^5 + /p^p(2a2<T)„^ • (<T2a2)^5, (1) 

/575 - + /pV«/3 • Cr75- (2) 

Here and are constants that parametrize the two-body interaction in the density 
and spin channels, respectively, the superscripts pp and ph refer to the particle- 
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Figure 1: Sum of polarization tensors contributing to the response function of a 
superfluid. The diagrams b, c, and d are specihc to superfluid systems and vanish in 
the unpaired state. 


particle and particle-hole channels, while the quantities (Tap are Pauli matrices in the 
spin space. 

The low-lying excitations of neutron matter can be obtained by expanding the 
relevant response functions in small parameters of the theory, for example q/kp, where 
kp is the Fermi wave number and q is the magnitude of the momentum transfer. If 
we are interested in radiation processes, i.e., in the time-like kinematical domain, the 
quantity qvp/u is small, where vp is the Fermi velocity and u is the energy transfer 
(typically of the order of the temperature). 

As there are four distinct propagators in the theory of superfluid systems (the 
so-called Nambu-Gor’kov structure), the response function is given by four different 
contributions, shown as Feynman diagrams in Fig.[T] We adopt the standard notation 
of the diagram technique for superfluid systems, following Ref. [10] . The shaded (full) 
vertices appearing in the polarization loops can be computed, in turn, by resummation 
of inhnite loop diagrams, which account for the modihcations due to the interactions 
in the medium (see for example Ref. i)- To leading order in a small-momentum 
expansion, the temporal (00) and spatial (jj) components of density (D) and spin 
{S) response functions may be written as 
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( 4 ) 


where 


g{v,u,q) 


+00 



—oo 


/( 6 -) - /( 6 +) 

- (e+ - e_y + ir] 


+ ^- 1 ~/(^-) ~/(^+) 
e+e_ - {e++ + ir] 


Here f{x) = {exp[(a; — /i)/T] -l- 1}”^ is the Fermi function, jj, being the chemical 
potential. We dehne e± = + and qv/2, where is the spectrum 

in the unpaired state, A is the pairing gap, v is the particle velocity, and q is the 
momentum transfer. In concentrating on the pair-breaking contribution, we need 
only keep the term proportional to 1 — /(e_) — /(e+). 
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Figure 2: Comparison between numerical (solid lines) and perturbative (dashed lines) 
results for the density response function 11^, in units of the density of states v [8]. 
The imaginary part of the response function is shown by a heavy (blue) line, the real 
part by a light (cyan) line. The energy transfer ca is in units of cjq- The temperature 
is 0.5 Tc and the pairing gap is A = 1.0 MeV. The ratio of momentum transfer 
and Fermi momentum is kept hxed at qjkp = 0.01, and the Fermi momentum is 
kp = 1.0 fm~^ (which corresponds to a density n = 0.22no in terms of the nuclear 
saturation density uq.) 


Figure E] illustrates the dependence on the transferred energy cu, of the real and 
imaginary parts of the density response function of neutron matter, for hxed three- 
momentum transfer. The zero-temperature gap is related to via = A(0)/1.76. 
The Landau parameter is set to = —0.5 for density perturbations, while the 
frequency and momentum transfer are normalized to the threshold frequency Uq = 
2A(T). The response function in the negative energy range can be obtained from 
the relations Ren'^(—ca) = Ren^(a;) and Imn^(—ca) = —Imn^(a;). There exists a 
threshold in the response function because an energy a; > cjq is needed to break a 
Cooper pair, i.e., the energy transfer must exceed the binding energy of the pair. For 
energies above the threshold (ca > cuq), a non-zero imaginary part implies that the 
collective excitations have hnite life-time. They are not perfect quasiparticles. 


3 Numerical results and collective modes 


The numerical results reported here were obtained in Ref. [8]. Once the response 
function is determined, we are in a position to write down the associated spectral 
function, which is dehned by R{u, q) = — 2Imn(a;, q) and can be cast as 


R{u,q) 


2{v^J-HmP{u,q) 

-ReP(a;,q)]%ImP(a;,qr)2’ 


( 5 ) 
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where P{u,q) = Q'^{u,q)/C{u,q) and 


q) = q)C{uj, q) - B{uj, q)V^{uj, q). 


The elementary loops are given by [mca 
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uPqv u-qv^ 
dA 


g{v,uj,q), 


(7) 


( 8 ) 


(9) 

( 10 ) 


where the operator V is equal to +1 for vertices which are even under the time reversal 
operation and —1 for vertices odd under this transformation, while u is the density 
of states. 

For small imaginary parts the spectral function can be approximated as 


R{uj,q) = 2TrZ{q)5{{fll) ^ - ReP{u, q)) + R,eg{uJ, q), (11) 

where Rreg{uj,q) is the regular (he. smooth) part of the spectral function and Z{q) 
is the wave-function renormalization [8]. The Dirac delta gives the spectrum of the 
quasiparticle excitations. In Eq. (ITT]) the interaction term indicates that one is 
dealing with density perturbations. 

The dispersion relations of the excitations are determined by solutions oj{q) of the 
equation 


1 - /SReP(u;, q) = 0. (12) 

Fig. |3] demonstrates its solutions in the case of density and spin excitations for the 
range of interaction values — 1 < <2. It turns out that if the interaction has a 

positive sign the modes are located in the domain cu/cuq > 1, where lmP{u),q) ^ 0. 
This is the case for spin excitations in the crust of a neutron star. Such modes are 
termed diffusive. The density modes are located in the domain uj/uo < 1, since 
the particle-hole interaction takes negative values in neutron-star crusts. The pair¬ 
breaking part of lmP{u,q) vanishes in this domain. These modes can be called 
excitonic in that they can be viewed as bound pairs of particles and holes (and 
are analogous in this respect to modes occurring in ordinary semiconductors). It 
must be noted that we cannot follow the high-momentum behavior of the modes 
because we are using perturbative response functions, which are no longer valid at 
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Figure 3: Dependence of the frequency of the modes on momentum transfer for 
density (left panel) and spin-density (right panel) perturbations [8]. The values of 
particle-hole interaction /ph are shown in each panel. The system parameters are 
hxed aX kp = 1 fm“^, A = 1.0 MeV, and T/Tc = 0.5. Heavy lines (blue online) 
correspond to undamped excitonic modes. Light lines (cyan online) correspond to 
diffusive damped modes. 


large momentum transfers q/kp ~ 0.3. We also observe that for negative /ph, the 
modes tend toward zero with increasing momentum transfer. The high-momentum 
behavior of the excitonic modes requires further study in order to realize a complete 
picture of their thermodynamics. In the following section we will concentrate on the 
properties of the spin diffusive modes. 


4 Thermal conductivity of spin diffusive modes in 
neutron-star crusts 


We conclude with an estimate of the thermal conductivity of spin diffusive modes, 
assuming that this property is dominated by the mode-mode scattering process. To 
treat the thermal transport in the presence of a thermal gradient VT, we write the 
kinetic equation for spin diffusive modes as 


de dfa 
dp&T 


vr = /[/], 


(13) 


where /o is the equilibrium distribution of the modes, which can be taken as a Boltz¬ 
mann distribution owing to the large value of Uq and /[/] is the collision integral. 
Assuming a small perturbation, we linearize the kinetic equation in the standard 
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fashion, i.e., we substitute 


/ = /o + if, Sf = ~X = JX- (14) 

in the collision integral. For this generic form of the perturbation, the collision integral 
has the form 

/[/i] = ^j w'fo2{x'i + X2 - XI - X2)dT[dT2dT'„ ( 15 ) 

where (iFi etc. denote the phase space volume measure and w' is the transition prob¬ 
ability. For thermal conductivity the perturbation is of the form 


^ = g.VT = vVTj2 ( 16 ) 

S=1 


where v = (de/dp), (3 = (4aT)“^, and the S^{x) are associated Laguerre polynomials. 
Truncating the expansion in polynomials by keeping only the leading s = 1 term, the 
thermal conductivity is given by k = 75/16aii where 


Oil — 


2^0 
T 


25/27rV2 


dv v'^e 


sin^ a 


da 

da 


da, 


(17) 


in which da/da is the mode-mode scattering cross section. It is expected that the 
contribution from terms of higher order in s is small. 

We can make a simple estimate of k by assuming a constant cross section given 
by da/dVt ~ 1/Ak\. We hnd 


K ~ 


2aT 


TT 


4' 


T 


(18) 


which exposes the temperature scaling of the thermal conductivity when the cross 
section is independent of T. 

In closing, we note the analogy of the spin diffusive modes with rotons in superfluid 
helium, which can provide a guide for further studies of scattering and transport 
involving these collective modes. 
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